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Abstract: Recent searches for first-generation leptoquarks and heavy right-handed 
Wr bosons have seen excesses in final states with electrons and jets. A bizarre 
property of these excesses is that they appear to violate lepton universality. With 
these results in mind, we study the phenomenology of supersymmetric models in 
which the Higgs arises as the sneutrino in an electron supermultiplet. Since the 
electron is singled out in this approach, one can naturally account for the lepton 
flavor structure of the excesses. In this work, we show that in such a framework, 
one can significantly alleviate the tension between the Standard Model and the data 
and yet evade current constraints from other searches. Lastly we point out that 
correlated excesses are expected to be seen in future multilepton searches. 
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1 Introduction 

The Standard Model (SM) of particle physics is among the most successful models 
ever devised, yet it leaves open several puzzles that should be resolved by a more 
complete description of nature. A well-motivated, broad class of models based on 
supersymmetry (SUSY) has the potential to resolve one or more of the outstand¬ 
ing puzzles of the SM, including the hierarchy problem, the nature of dark matter, 
the mechanism of baryogenesis, and the running of gauge couplings to a unified 
value. From a phenomenological point of view, however, there are several issues with 
models based on SUSY. In particular, the naive implementation of natural R-parity 
conserving MSSM requires a light spectrum of color-charged particles to which the 
experiments at the Large Hadron Collider (LHC) should have sensitivity, yet no hints 
of SUSY have been seen in the “standard candle” channels with Missing Transverse 
Energy [1], Furthermore, a Higgs boson with mass 125 GeV is not generically rec¬ 
onciled with a natural spectrum of superpartners [2]. Both of these tensions hint 
at the possibility that, if natural SUSY describes our universe, then it may have an 
alternative structure. 

The lack of observation at colliders has lead to the introduction of many vari¬ 
ations of supersymmetry such as R-parity violating (RPV) [3-11] and A-symmetric 
supersymmetry [12-22, 22-33]. Constraints on SUSY, even in the context of RPV 
models, are already quite stringent [34-36]. These constraints are somewhat less re¬ 
strictive in models with R-symmetric models due to the requirement of Dirac gaug- 
inos [37]. In particular, this prevents same-sign lepton signatures that would be 
smoking gun indicators of physics Beyond the SM (BSM). An additional intriguing 
feature of such models is that they allow for the Higgs field to be identified with 
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the superpartner of a left-handed electron [22, 38, 39 ] 1 . In this unique framework, 
traditional LLE C and LQD C RPV effects are present but necessarily suppressed by 
the smallness of the Yukawa couplings. However, RPV effects appear due to a mixing 
between the electron doublet and the gauginos (such mixing has been previously used 
to put constraints of possible sneutrino VEVs [40]). Since the electron is singled out 
as the Higgs partner, such models have non-standard lepton flavor structure leading 
in general to an abundance of electrons in the final state. Furthermore as we will 
show, the requirement of nearly massless neutrinos requires the introduction of an 
R-symmetry. 

The CMS experiment has recently seen hints of potential BSM physics at the 
~ 2.5a level in three separate searches that appear to single out the first generation 
of leptons. Two of these analyses were optimized to look for pair production of 
leptoquarks. In one case, the leptoquarks decay to an eejj final state, while in the 
other they decay to an eujj final state [41] . Both showed excesses hinting at a roughly 
650 GeV leptoquark, at the 2.4a and 2.6a levels respectively. The excesses are not 
consistent with the only decay modes of the leptoquarks being ej and vj [41, 42], 
The third search was optimized for a Wr decaying to an eejj final state and saw a 
2.8a local excess for a resonance near 2.1 TeV [43]. However, the distributions of the 
excess do not appear to be consistent with those of a naive Wr [43] (though see [44] 
for a more general discussion on this possibility). Its important to note that the 
leptoquark searches did not see an excess in its high leptoquark mass bins. While 
not emphasized in earlier work, this puts serious limitations on new BSM signals 
attempting to explain the excess. No excesses were observed in the corresponding 
channels with muons [43, 45]. 

Several models have been constructed in order to explain this excess. Many 
of these models are supersymmetric in nature [46-51] (see [42, 44, 52-63] for non- 
supersymmetric explanations). The vast majority of them do not attempt to explain 
the puzzling lepton flavor structure of the observed excesses, but merely choose 
certain couplings to be larger then others. Standard tools for suppressing flavor- 
violating processes such as minimal flavor violation (MFV) [64] cannot explain a 
different coupling for the first and second generations. In MFV, such non-universal 
terms in the Lagrangian are suppressed by m M /m T . Furthermore, due to the presence 
of a heavy resonance, these models often predict an excess in the searches for higher 
mass leptoquarks, which has not been observed in the data. 

In this paper, we investigate the possibility that supersymmetric models with 
the Higgs as a sneutrino could explain the excesses seen by CMS. The lepton flavor 
structure is naturally obtained within the context of such models. The complex SUSY 
spectrum yields a rich variety of decay modes, suppressing the number of events 

1 In general this can be any lepton, but as we will discuss in section 2, the electron is the most 
natural choice. 
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seen in individual channels and allowing such models to evade many constraints. 
Overall, this class of models provides a good fit for the current data, while making 
several new and testable predictions for the upcoming run of the LHC. The role 
of the leptoquarks in the model is played by a left-handed first-generation squark 
with R-parity violating decays, while the heavier ~ 2 TeV resonance is explained by 
gluino-squark production. The masses that give the best fit are an up squark mass of 
810 GeV and a gluino mass of 1790 GeV. In addition to accounting for the excesses 
observed by CMS, this model addresses the lack of an excess when the set of cuts is 
optimized for higher mass leptoquarks. The model considered in this paper addresses 
this potential issue by softening the “leptoquark” spectrum with additional jets, as 
proposed in [60]. 

The remainder of this paper is structured as follows. In Section 2, we review the 
minimal model with the Higgs as a sneutrino. We determine a set of parameters of 
this model that provide a good fit to the current CMS data in Section 3. We then 
conclude discussing current bounds on the model and provide additional predictions. 


2 Model with Higgs as a slepton 


2.1 Overview 

To illustrate the main ideas behind the Higgs-as-slepton model [38], we begin by 
attempting to construct a supersymmetric Standard Model that is more minimal 
than the MSSM. One can identify the SM Higgs doublet H with a slepton doublet 
L a , since they are both in the same gauge representation (1,2)_!/2- The model 
then requires two fewer doublet chiral superficlds than the MSSM. However, a major 
issue arises from the fact that the Kahler potential generates electroweak-scale Dirac 
masses between the partner leptons L a = (u a , l~) and the Winos and Binos: 




, 9 V H 5 . , 

H-—+ h.c. 


( 2 . 1 ) 


This leads to neutrino masses that are too large. 

One way around this difficulty is to first impose a U (1)# symmetry, with R-charge 
zero for the slepton doublet L a and —1 for the partner lepton doublet L a . The U(1)r 
symmetry remains unbroken when L a acquires a VEV, and can still forbid Majorana 
masses for all C/(l)#-charged neutralinos. By introducing adjoint chiral superfields 
$ and SUSY-breaking Dirac gaugino masses, one of the neutralino mass eigenstates 
becomes massless. This massless neutralino is mainly comprised of v a and can be 
identified with the “physical” neutrino. 

We now present the details of the model. Table 1 lists all the superfields and their 
gauge and U(1)r representations. With the CMS excesses in mind, we have chosen 
the first-generation leptons to partner the Higgs. This will give rise to experimental 
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Table 1: Superfields and their gauge and U(1)r representations in the Higgs-as- 
slepton model. 


signatures specific to the electron, without the need to tweak any lepton couplings. 
B and L are arbitrary parameters that determine the U( 1)# representations of the 
quark and the 2nd- and 3rd-generation lepton superfields. 

The most general superpotential consistent with the symmetries (assuming B ^ 
1/3 and L ^ 1) is 

3 2 

W = !Jd,ijHQ l D c - + VejHLiE?. (2.2) 

i,j=1 Z—1 


We have chosen to work in the mass basis of the charged leptons. The superpotential 
does not generate up-type quark masses due to the absence of an up-type Higgs 
superfield H u . The same is true for the electron mass, since the required term 
HHE 3 is identically zero. Both can be generated by SUSY-breaking Kahler terms 
of the form [38] 


J d 2 dd 2 e 


Xt HiQiU? 
M A 


(2.3) 




(2.4) 


that are suppressed by a A cutoff scale. This also provides a natural explanation for 
the smallness of the electron mass, hence further motivating our decision to partner 
the first-generation leptons with the Higgs. 

The U (l)ij symmetry forbids mixing between left-handed and right-handed squarks, 
so the squark phenomenology differs from that of the MSSM [25]. This also simpli¬ 
fies our subsequent analysis of squark production and decay since the squark mass 
eigenstates are then either left- or right-handed. 

We note that the terms in the superpotential can also be interpreted as RPV 
terms of the form L^QiD^ and L^LiEj. Therefore, experimental bounds on RPV 
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coefficients [65] can be applied to the superpotential Yukawas y ( uj and y e ,ij, which 
are in turn determined by the SM fermion masses and mixings. We fold that these 
bounds are satisfied by the model for the choices of squark masses to be used in later 
sections. 

While we assume the model described above in this work, our results are largely 
independent of the detailed mechanism giving the up-type quark and electron masses. 
Alternative models which introduce additional chiral superfields are also possible [22, 
39] and can also produce similar signatures. 

2.2 Chargino and neutralino mass matrices and mixing 

The chargino and neutralino Dirac mass matrices are given by 
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We have neglected the masses from A-suppressed SUSY-breaking terms such as elec¬ 
tron masses, since they are much smaller than the present terms and hence not 
expected to play an important role. To order e = g%/(2M^) = mw/M the 
chargino 4 component mass eigenstates are: 
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( 2 . 6 ) 


with mass eigenvalues m- = 0 and m- = m- — AAp. The mass eigenstates for 

Xp X2 X3 vv 

the neutralinos are 2 : 


Xi = 


0 


Y 0_fC + e ^l yO _ 
X'2 1 I ’ X3 


i_ _ g’ 

US g tf- We 


B c 


(2.7) 


with mass eigenvalues m x 0 = 0, m x 0 = and m x 0 = Mg. 

X\ can be identified with the physical electron, and Xi with the “physical” 
electron neutrino, before PMNS mixing. We note that the gauge couplings of the 
physical gauginos and first-generation leptons to W ± and Z are affected by the 0(e) 
mixing. One consequence is that the eeZ coupling is modified, hence violating lepton 
flavor universality. This allows us to place a lower bound of ~ 2 TeV on the Dirac 


2 We have assumed here that 


M- - M 2 - 

1V± W B 


M 2 — M 2 

W B 


< 


rriw- In the converse case where 
the actual heavy neutralino eigenstates are linear superpositions of \2 an< l X3 above, with 
mixings given by the Weinberg angle 9\y- Nonetheless, this does not affect any of our subsequent 
results on the partial widths. 
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Figure 1: Mixing-induced decay channels in which a supersymmetric particle qp 
decays completely to SM particles. 


chargino mass [38]. Another consequence is that the modified gauge couplings 
mix the physical gauginos and leptons, thus providing a channel for the gauginos 
to decay completely to SM particles, e.g. X 2 ~► Xi W + . Should the squarks be 
lighter than the gauginos, which we assume in the rest of this work, virtual cascades 
such as cLl — y d X2 ~t d x± IF - may also become important decay channels for the 
first-generation squarks, as we will see below. 

2.3 First-generation left-handed squark decays 

In MSSM with RPV, supersymmetric particles can decay completely to SM particles 
through channels generated by RPV superpotential and soft SUSY-breaking terms. 
While this is also true for the Higgs-as-slcpton model, there are new decay channels 
due to the mixing of physical gauginos and leptons by the modified gauge couplings. 
A typical diagram for the new channel is shown in Fig. 1. The new channels are 
especially important for first-generation squarks compared to the standard RPV 
channels, due to the smallness of the Yukawas in the latter [38]. The approximate 
partial widths of these channels for first-generation LH squarks are shown in Table 2. 
Fig. 2 compares the partial widths of the mixing-induced and standard RPV channels 
for di decay, from which we see that the former is dominant except for very large 
values of 

Supersymmetric particles (and the Higgs) can also decay into SM particles + 
the gravitino, which is the lightest supersymmetric particle (LSP) in the model. The 
decay occurs via goldstino interaction terms fixed by supersymmetry, with partial 
widths that typically scale as ml p / (m 3 / 2 Mpi) 2 , where m sp is the sparticle mass, 1713/2 
the gravitino mass and Mp\ the Planck scale [38]. However, as long as the gravitino 
mass is not too small (m 3 / 2 leV), these decays are expected to be sub-dominant 
and can hence be neglected. For the rest of this work, we assume all first-generation 
squarks to decay via the mixing-induced decay channels. 
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Decay channel 
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Table 2: Partial widths for the mixing-induced decay channels. Here xi and Xi 
refer to the physical electron and electron neutrino. Yq is the hypercharge of the LH 
quark doublet. The decay channels have been arranged such that the approximate 
isospin symmetry from the Goldstone boson equivalence theorem is obvious. 



Figure 2: Partial widths of cR for mixing-induced and standard RPV decay chan¬ 
nels, assuming m d - = 810 GeV and Mg = M^. The mixing-induced channel domi¬ 
nates over the range of M^ considered. 


3 Simulation and Results 

In this section, we estimate the contribution of the above model to the CMS lep- 
toquark and Wr searches. The spectrum and production channels of interest are 
depicted in figures 3 and 4. 

The model predictions are calculated at tree level using Madgraph [66], Pythia 
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Figure 3: 


The spectrum of our benchmark point. All other fields are decoupled. 
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Figure 4: Sample production mechanisms for disquark and single gluino production 
channels. Squarks decay through the 3 body decay shown in figure 1. 
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Figure 5: The m ee jj prediction for our model after applying cuts used in the Wr 
search. The background and relevant cuts were taken from [43]. 


6.4 [67] for showering and hadronization, and PGS [68] for detector simulation. The 
model files were created using Feynrules [69]. To estimate the next-to-leading order 
(NLO effects we scaled the cross-sections by their corresponding K-factors calculated 
using Prospino 2.1 [70]. While Prospino was designed for the MSSM we do not expect 
significant deviations in the calculations of K-factors. 

The Wr search distribution is shown in figure 5. We reproduce the invariant 
mass distribution of the two leading electrons and two leading jets. We also applied 
all the relevant cuts detailed by CMS in Ref. [43], the most restrictive requiring the 
invariant mass of the electrons be greater then 200GeV. 

































(a) eejj search (b) evjj search 


Figure 6: Bin-by-bin background-subtracted events for the LQ searches. Each bin 
count is a subset of the previous bin and hence the bins are highly correlated. The 
model shows some tension with the data at high LQ mass cuts. 


The single gluino production dominates the high mass peak, while the disquark 
channel contributes broadly to the bins between 1 — 2 TeV. The broad feature is 
a consequence of a many-body structure of the decay which, and as pointed out 
in [60] , is useful to evade bounds by the CMS leptoquark search without introducing 
multiple decay channels. We emphasize that in our model we satisfy both properties 
of the signal. Firstly, no signal is found in corresponding muon channels as only the 
electron doublet mixes with the other neutralinos and charginos in this framework. 
Secondly, the events are dominated by opposite-sign electrons. This is guaranteed by 
the imposed R symmetry for which an electron and positron have opposite charges. 

Next we reproduce the leptoquark (LQ) searches in this framework. In the 
LQ search a sequence of more stringent cuts are applied, optimized for different 
mass leptoquarks. In the eejj channel, the main discriminating variables are St 
(the scalar sum of px of two leading electrons and jets), m ee (invariant mass of the 
two electrons), and m™ n (the minimum of the electron-jet invariant mass of the 
four possible combinations for eejj). In the evjj channel, the main discriminating 
variables are St, E™ lss , and rn ej . Typically models that predict large rri eerj (in order 
to explain the Wr excess) will also produce large St (and rn 7 ™' n unless they arise from 
a very light LQ). In general, this leads to expected excess in the heavy LQ mass cut 
range. Thus it is important to check the predictions of any model attempting to 
explain the CMS excesses in these high mass bins. 

The corresponding cuts for each LQ mass can be found in [41] (see tables 2 
and 3). Here we plot the difference between the data and the SM background as a 
function of LQ mass cut. The results are shown in Fig. 6. Each bin is a fraction 
of the events in the lower LQ mass cut bin and thus the bins are highly correlated. 
We see moderate agreement of our signal with the observed counts. We are able 
to explain the excess in the ~ 650GeV region, but see small excess in the higher 
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(a) eejj search (figure 5 in [41]) 


(b) evjj search (figure 8 in [41]) 


Figure 7: The CMS leptoquark search plots. 


mass cuts for eejj. The excess in the high mass range is a general characteristic of 
trying to explain both the Wr and LQ searches. Note that the excess is 0( 5) events 
instead of (9(10) which were found in the Wr search. This is a consequence of the 
large number of jets increasing the effectiveness of the Sr cut. 

To further check the kinematic properties of the model we compare our m 7 ™' n and 
m e j distributions at the 650 GeV mass cut point. The results for both searches are 
shown in figure 7. In both the eejj and evjj channel we see good agreement between 
the model and experiment. The broad feature of the plots is again a consequence 
of the many jet signal and is necessary to get the right kinematic spread in the LQ 
invariant mass distributions. 

This framework has two characteristic features - many electrons in the final state 
and many jets. Due to their limited background, we expect the most stringent bounds 
on our model arise from multilepton searches [71, 72], The model produces more than 
2 leptons if each squark decays into an electron and additional leptons arise from 
vector boson decays. We now roughly estimate the number of expected events in the 
multilepton searches. The NLO cross section for squark-squark and squark-gluino 
production at our mass point is 5.7 fb. At C ~ 20 fb _1 this corresponds to about 
115 events. The probability of both squarks producing electrons (as opposed to 
neutrinos) is about 1/4. Furthermore, the probability of at least one of the vector 
bosons decaying leptonically is between 11 and 40% depending on whether there is 
a WW, WZ, or ZZ is in the final state. This suggests 5 — 10 events with 3 or more 
leptons. However, these events don’t contain any genuine E™ lss or 6-tagged jets, 
both of which are powerful discrimating variables in such searches. This makes the 
signal hard to detect, even in a multilepton search. Thus we conclude the model is 
safe from current multilepton bounds, though we expect sensitivity with more data 
at higher energies. 
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4 Discussion and Conclusions 


In this paper, we have explored the phenomenology of a class of SUSY models in 
which the Higgs is a sneutrino. Such models could account for excesses seen in the 
CMS experiment, while accounting for the observed kinematics and flavor structure 
in a natural way. 

As with most SUSY models, several correlated observables are expected. While 
the detailed spectrum and branching fractions are model-dependent, these models 
have a few generic predictions. Most reliably, there should be correlated excesses in 
multi-lepton searches. Since the decay of hadronic sparticles necessarily proceeds via 
clectroweakinos, the decays will generally feature leptons, possibly in large numbers 
and with a preference for electrons. These excesses would come with some missing 
energy from neutrinos, but decays without neutrinos are certainly possible. The 
lepton structure of these excesses would again be striking, featuring more electrons 
than muons or taus. The scales of < 1600 GeV from qq*, < 2400 GeV from qg, 
and < 3600 GeV from gg would also feature in the total invariant (transverse) mass 
distribution. 

The remaining signals are highly dependent on the more weakly coupled or heav¬ 
ier elements of the spectrum. The constraints on slcptons and electroweakinos remain 
weak after Run 1 of the LHC, but searches for signatures of new electroweak states 
are a vital part of Run 2 that can only be fully exploited at high luminosity. Such 
particles with mass (9(100 GeV) could be in the spectrum and would decay primarily 
to clecrtoweak bosons, electrons, and neutrinos. 

The first run of the LHC has seen a remarkable confirmation of the SM with few 
searches finding excesses beyond the 2 o level. On the other hand, several searches 
that have seen excesses indicate similar final states with electrons and jets, as well 
as large energy scales of ~ 650 GeV and ~ 2 TeV. If such excesses are the first hints 
of a new state beyond the SM, then Run 2 will bring striking and nearly immediate 
discoveries, as the sensitivity to physics at ~ 2 TeV is vastly superor to that in the 
first run. 
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